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Abstract 

The chiral symmetry of QCD is studied at finite temperature and chemical 
potential using the Schwinger- Dyson equation in the improved ladder approx- 
imation. We calculate three order parameters; the vacuum expectation value 
of the quark bilinear operator, the pion decay constant and the quark mass 
gap. We have a second order phase transition at the temperature T c = 169 
MeV along the zero chemical potential line, and a first order phase transition 
at the chemical potential \x c = 598 MeV along the zero temperature line. We 
also calculate the critical exponents of the three order parameters. 
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1 Introduction 



The chiral symmetry in QCD is dynamically broken at zero temperature. This feature 
is confirmed by the fact that the pion is the Nambu-Goldstone boson accompanied with 
this symmetry breaking. On the other hand, it is shown that spontaneously broken 
symmetries restore at sufficiently high temperature (and/or chemical potential) in some 
simple models. [0] Then, the same restoration is also expected to hold for the dynamically 
broken chiral symmetry in QCD. This phenomena is widely believed to be seen in heavy- 
ion collisions, the early universe and the neutron stars. 

There are various attempts to study the phase diagram and critical behavior. In order 
to do them we need any non-perturbative treatments such as e or 1/N expansion, lattice 
simulation, the Schwinger- Dyson equation and so on. Based on universality arguments 
it is expected that critical phenomena of finite temperature QCD is described by a three 
dimensional linear a model with the same global symmetry, in which the e expansion is 
used.[U Lattice simulations are powerful tools to study the QCD at finite temperature. 0, 
[|, |5|, H| The Schwinger- Dyson equation in the improved ladder approximation is solved 



with further approximations and give the dynamical symmetry restoration. M, 0,0, 10, 11 



Nambu-Jona-Lasinio models, as phenomenological models of QCD, provide us with useful 



pictures about the dynamical chiral symmetry breaking and its restoration. 0, [13], |Tj] In 
these three approaches it is indicated that there is a second order phase transition at 
(T, /i) = (T c , 0) and a first order one at (T, /i) = (0, /i c ) in the case of the two massless 
flavors. The phase transition points are found to be of order T c ~ 200, /i c ~ 400 MeV. 

In this paper we use the Schwinger-Dyson equation in the improved ladder approx- 
imation. The advantages of this approach are that it is a convenient tool to study the 
nature of the chiral symmetry, and we easily introduce fermions coupling to the gluon in 
the chiral limit at finite temperature and chemical potential. Further, we have no degrees 
of freedom (parameters) to fit the physical observables, and we obtain a definite answer. 

In previous attempts further approximations are introduced in addition to the ladder 
approximation. Some non-perturbative approximation can violate the chiral symmetry, 
and reliable results cannot be obtained. We would like to obtain the results keeping 
the chiral symmetry within the framework of the Schwinger-Dyson equation only in the 
improved ladder approximation. 

We calculate the values of the three order parameters; the quark mass gap, the vacuum 
expectation value of the quark bilinear operator (ipip) and the pion decay constant. We 
have a second order phase transition at T c = 169 MeV along the fi = line and a first 
order phase transition at /i c = 598 MeV along the T = line. The critical exponents of 
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the above three order parameters are extracted at (T, /i) = (T c ,0), which shows that our 
formulation is different from mean field theories. 

This paper is organized as follows. In section 0, we show basic ingredients to study 
the chiral symmetry restoration using the Schwinger-Dyson equation in the improved 
ladder approximation. The expressions of the three order parameters are given in terms 
of the quark mass function. The Pagels-Stokar formula is used to calculate the pion decay 
constant. In section |3] we give our numerical results. The Schwinger-Dyson equation is 
numerically solved using a iteration method. We determine the positions and the orders 
of phase transitions. We also extract the critical exponents. Summary and discussion are 
found in section |j. 

2 Schwinger-Dyson Equation at Finite Temperature 
and Chemical Potential 

The restoration of spontaneously broken symmetry occurs at finite temperature and chem- 
ical potential. The phenomena is described in terms of the imaginary time formalism 
in gauge theories. [15|, [16| 

In this section we show basic ingredients to solve the Schwinger-Dyson equation at 
finite temperature and chemical potential. Then, we study the dynamical chiral symmetry 
breaking and its restoration. In this paper all dimensionful quantities are rescaled by using 
the Aqcd, otherwise stated. 

We consider the QCD with massless u and d quarks, and then there is a SU(2)l x 
SU(2)r chiral symmetry. There are several probes to investigate the chiral symmetry 
such as the quark mass gap, the vacuum expectation value of the quark bilinear operator 
(VEV) and the decay constant of the pion. Those are evaluated in terms of the quark mass 
function The mass function is determined by the Schwinger-Dyson equation. We 

use the improved ladder approximation to solve this equation. We work with three flavor 
/3-function for the running coupling, since the s quark also contributes to the running of 
the coupling in the concerned energy range. 

At zero temperature this approximation provides us with a convincing result of the dy- 
namical chiral symmetry breaking and good values of the lowest-lying meson masses. So, 
we expect that the approximation gives good result even in the case of finite temperature. 

To take the effect of finite temperature into account, we work in the imaginary time 



formalism |15|, [L6|. Let us start with writing down the Schwinger-Dyson equation as in 
Fig. [1|. The diagram is exactly the same as that in the zero temperature QCD, since the 
difference between the usual (zero temperature) field theory and finite temperature field 
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Figure 1: The Feynman diagram of the Schwinger-Dyson equation. We 
have the same diagram for the equation at finite temperature and chemical 
potential. 



theory stems from the boundary effect of time only. The time components of quark and 
gluon momenta become discrete. Since quark fields have anti-periodic boundary condition 
in the imaginary time direction, we have 

p° = 2niT (n+- 
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k° = 2mT(m+-) , (2.1) 

where n,m e Z. When the chemical potential \x is introduced, the time component of p 
in the quark propagator Sp(p) is modified as 

Po -> Po - At . (2.2) 
The momentum integration is modified to the summation: 

Then, modifying the Schwinger-Dyson equation at zero temperature according to Eqs. fl2.1|) , 
and (|2.3|), the Schwinger-Dyson equation at finite temperature is given as 



^-S F {p)- l =T E ^C 2 g 2 (p,k)K^(p-k)rS F (k)Y , (2.4) 

m=-oo J \ A7Y ) 

where C 2 is the second Casimir invariant and —K^ v is the gluon tree propagator in the 
Landau gauge: 

K,u(l) = - l -f) ■ (2.5) 

The quantity g 2 (p,k) is a one- loop running coupling depending on the momenta p and 
k. In order to describe a property of QCD we use the following form for the running 



3 



coupling: JL7 



9 2 (p, k) 



0o 



x < 



1 

7 
l 

l 



(t F -t c ) 2 -{t-tcf 



+ 



2t 2 F {tp 
{t F ~ t c ) 



tc) 



2t% 



if tp < t 

if tc < t < tp 

if t < t c 



(2.6) 



where t = ln(p 2 + k 2 ), t c = -2.0, f3 = (llN c -2N f )/(48n 2 ) is the coefficient of one-loop (3- 
function and tp is a parameter needed to regularize the divergence of the running coupling 
at the QCD scale Aqcd- We call tp the infrared regularization parameter. N c = 3 and 



N 



f 



3 are the number of colors and flavors, respectively. The running coupling g 2 is 
smoothly interpolated between the ordinary one- loop running coupling form at t > tp and 
a constant value in the low energy region. As will be shown, the results do not depend on 
this particular infrared regularization. By virtue of the running effect of the coupling the 
resultant mass function, which is determined by Eq. (|2.4j), reproduces the exact behavior 



in the high energy region. [|T^] Notice that this property is needed to preserve the chiral 
symmetry. [ |19| 

The quark propagator is expanded by three SO (3) invariant amplitudes as 

1 



Sf(p) E(p) + ( f i + B(p))^-A(p^- 

At the vanishing temperature and chemical potential limits (T, /j, 
Landau gauge allows us to obtain 



(2.7) 
0) the choice of 



A(p) 
B{p) 



1 , 
. 



(21 



Although we are studying in finite temperature and chemical potential, we assume the 
relations (|2.8| ) for simplicity. We expect that the relation (|2.8|) is not changed so much in 
the case of low temperature and chemical potential. As shown later, the phase transition 
line of the chiral symmetry restoration lies in that region; T c , \x c ^Aqcd- Then, the result 
should not change qualitatively. 

Now, substituting Eq. (|2.7|) into Eq. (|2.4|) under the condition Q2.8|) , we obtain 



53 / ydyK nm {x,y) 



m=— oo 



27rT(m+±)+m) + y 2 + Y> m (y) 2 



(2.9) 



where x = \p\ and y = \k\ and 



v , . 3TC 2 g 2 (p,k) (A<K 2 T 2 {n-m) 2 + {x + y) 2 " 
K nm (x,y) = — In 



8w 2 xy 



\47r 2 T 2 (n 



m) 



x 



yf 



(2.10) 
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Notice that we have the SO (3) rotational invariance but not 5*0(3, 1). The mass function 
S(p) is a function of p° and p 2 , and is rewritten as T, n (x). In the presence of the chemical 
potential we easily find from Eq. (|2.9| ) that the mass function takes a complex value 
satisfying the relation 

£_ n (.x)* = £n-i(x) for n=l,2,.... (2.11) 

We solve the Schwinger-Dyson equation (|2.9|) numerically by an iteration (relaxation) 
method. The momentum valuables x, y are discretized to be 

x -> x w = exp (A Jfl + (A yv - A//;) ~~ , (2.12) 

V ^SD — 1 / 

where n = 1, 2, ■ • • , and similarly for y. We divide In a; and hiy into Afe^ points. The 
quantity A = exp Kuv defines the ultraviolet cutoff for the space component of momenta. 
Therefore an 5*0(3) symmetric cutoff A is introduced, which is needed for numerical 
calculation. The momentum region of the time component is properly truncated so that 
the support of the mass function is covered well, as well as that of the space component. 
We have integrable singularities at (n, x) = (m, y) stemming from the tree level gluon 
pole, but these should be regularized in the numerical calculation. In order to avoid this 
singularity we apply a two-point splitting prescription 

K nm (x,y) -> ^(K nm (x,y + ) + K nm (x,y^)) , (2.13) 

with y± = y exp(±(A uv — A IR )/ (4N S d))- The validity of this prescription is checked by 
using the conventional (zero temperature) Schwinger-Dyson equation. 

After obtaining the mass function, we immediately evaluate the pion decay constant 
by using the Pagels-Stokar formula at finite temperature: 

x dEr,(x) 



/.(7>) 2 = ^F£/ * 2 dx- ^ 3 dx J (214) 



7T 2 *r* Jo 



We also calculate the VEV 



{27iT{n + \) + i>) +x 2 + E r 



fx) 2 



Ha = (dd) A = [ A * 2 dx- ^$ . (2.15) 

(2nT(n + ±) + ifA + x 2 + £„(x) 2 



vr 2 — Jo 



The VEV is renormalized at 1 GeV via 

11AT -27V 



/ 

- /InflGeVU ^ - 

(^)lGeV = \ nA J W>A , (2-16) 



where ip = u,d. 
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3 Numerical Results 



In this section we solve the Schwinger-Dyson equation numerically by an iteration method. 
We start with an initial form of the mass function and input it in the R.H.S. of Eq.( |2.9| ). 
Performing the integration of x and the summation of n, we have an updated form of the 
mass function, which is taken as a more suitable trial form. After sufficient iterations the 
functional form converges giving the true solution with enough accuracy. The convergence 
of the solution is very rapid off the phase transition regions. A typical form of the mass 
function is shown in Fig. ^| at T = 90, fi = MeV. The mass function dumps so fast in 




Figure 2: A typical form of the mass function. We put T = 90, /i = MeV. 
The integer n specifies the time component of the momentum as in Eq. ( j2.1|) 
and x is given in Eq. ( |2.12j ). 

the = —ipo direction that the values of the mass function at n ^ 4 are much smaller 
than that at n = (S n (x) ~ 10~ 2 x S (x)). This implies a dimensional reduction at 
sufficiently high temperature; i.e., four dimensional theory at finite temperature belongs 
to the same universality class as that of a three dimensional one with the same symmetry. 

We use the VEV and the pion decay constant as order parameters of the chiral sym- 
metry. The dynamical mass of the quark itself is also an order parameter. 

We determine the value of Aqcd from the experimental result f n = 93 MeV at T = 
\i = 0, and we have Aqcd = 592 MeV using Eq. ( p.l4j ) at tp = 0.5. This is our result 
of Aqcd obtained by using the one-loop /5-function. In this paper we put the infrared 
regularization parameter tp = 0.5 and show later that the physical observables as well as 
the phase transition line do not depend on tp. 
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3.1 Zero chemical potential case 

First, we study the phase transition along the /i = line. The phase transition point is 
defined so that three order parameters of the mass gap S n=0 (x = 0), the VEV (ipip) and 
the pion decay constant f„ vanish. We show the temperature dependences of these order 
parameters in Fig. ^| with two massless flavors. The SU(2)l x SU(2)r chiral symmetry 
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Figure 3: The functional forms of the order parameters (the mass function, 

the VEV and the pion decay constant) for T along the /i = line. 



restores at T = T c = 169 MeV. We have a second order phase transition. We have 
the same result as that of lattice simulations || § with two flavors in which the phase 
transition is second order at T c ~ 200 MeV. 

The ladder approximation, used here, gives no flavor dependence, since the depen- 
dence essentially comes through only the running effect of the coupling. Whereas the 
flavor dependence is suggested by the universality arguments^ and it is confirmed by 
lattice simulations II |5], ||, where we have a second order phase transition at Nf = 2 and 
first order ones at Nf > 3. The Nambu-Jona-Lasinio models pj, p~3| , I4f imply that the 
inclusion of the effect of U(1)a anomaly, so called the instanton effect, allows us to obtain 
the same flavor dependence as that in the lattice simulations. 

The important point in studying the chiral symmetry is that the approximation used 
should preserve the symmetry. Fortunately, the ladder approximation itself is consistent 
with the chiral symmetry. [|I9|1 However, many investigations violate the chiral symme- 
try, where further approximations are used in addition to the ladder M, |TU|. In order 
to reserve the symmetry the high energy behavior of the quark mass function must be 
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consistent with the result of the operator product expansion (QPE):|19| 

9C 2 



^ {fid n nx yx Nc -2N f 
X 



pi + p 2 ~ 00 



(3.1) 



Here we should notice that even in the finite temperature case the high energy behavior of 
the mass function is the same as that of zero temperature case, since the temperature effect 
is suppressed in the high energy region. In Refs. J7|, [TIJ their further approximation 
(X = const.) does not satisfy Eq. ( |3.1| ), on the other hand in Ref. JTT] they adopt an 
ansatz consistent with Eq. (|3.1| ) up to the logarithmic correction. While our formalism 



exactly reproduce the OPE result (|3TT 



We check the dependence of the order parameters on the infrared regularization pa- 
rameter tp. The physical observables, (ipip)iGeV and f n {T), should not depend on the 
parameter tp. We confirm this requirement. The dependence of the VEV and the f-^iT) 
are shown in Figs. f| and |5|. The values of the VEV and the /„- change, at worst, by 5% 
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Figure 4: The t F dependence of the VEV. It changes by 5% at worst against 

t F . 



and 8% against t F = 0.2 ~ 0.6, respectively. These values of tp = 0.2 ~ 0.6 correspond 
to those of the running coupling g 2 (p = k = 0) = 570 ~ 92.6. Moreover, the phase transi- 
tion point is fairly stable against tp. We conclude that there is no t F dependence of the 
physical observables and the position of the phase transition point. 

We also show the dependence of the mass function in Fig. ||. We have a strong t F 
dependence. Since the mass function is not a physical observable it may depend on the 
regularization parameter tp. However the phase transition point determined using the 
mass function is fixed. 
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Let us examine the critical behavior of the system. Since we have a second order phase 
transition at (T, //) = (T c , 0), the three order parameters behave near the phase transition 
point as 



W)lGeV ~ 
£ (0) ~ 



T\/3 
T 



1 - 



c 

7; 



(3.2) 



where T < T c . The critical exponents (3, v and /3' are numerically extracted by using the 
X 2 fitting. The order parameters O (= ( , 0V')iGeV) Sq(0), f^iT)) are fitted by the linear 
functional form lnO(T) = A + 7ln(l — T/T c ), and the fitting parameters are A and 7 
(= j3, u, (3 1 ) with T c = 169 MeV. We have good fittings, which is shown in Fig. [7|. The 
result is 

P = 0.171 , v = 0.497 , p' = 0.507 . (3.3) 

These values are different from those in mean field theories. If we consider mean field 
theories, we would have the relation (3 = v since £ ~ (ipip)- 
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Figure 7: The \ 2 fittings for extracting the critical exponents of the VEV, 

the mass gap and the pion decay constant. We draw the best fitted lines of 
the form A + 7ln(l - T/T c ). 
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3.2 Zero temperature case 

Next, we study the phase transition along the T = line. As is seen from Eq. ( |2.9| ), the 
mass function has an imaginary part for la ^ 0. The functional forms for fi are shown of 
the VEV, the pion decay constant f n (n), the real and the imaginary parts of the mass 
gap in Figs. || and |^. The SU (2) L x SU(2) R chiral symmetry restores at la — // c = 598 
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Figure 8: The functional forms of the VEV and the pion decay constant for 

(a along T = line. 
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MeV. We have a strong first order phase transition. Here, we check that the infrared 
regularization parameter tp does not affects the physical observables and the nature of 
the phase transition. 

Let us compare our result with those of other approaches. There is no lattice simulation 
at so large chemical potential that we can directly see the phase transition. A phase 
transition is suggested by extrapolating the result of lattice simulation around small /i.Q 
On the other hand, there are many attempts using Schwinger-Dyson equations f7|, || [TTJ 
and NJL modelsjT^, Tl|. The previous attempts^, [TTJ] in the ladder approximation give 



a first order phase transition. The NJL[TB|, [TJJ models with the instanton effect also give 
the same at Nf = 2, 3. Therefore, our result confirms theirs. Our advantage is that we 
have no parameter which modifies the physical result. 

3.3 the phase diagram 

Finally, we study the phase diagram of the chiral symmetry restoration. Near the /i = 
line we have second order phase transitions and near the T = line we have first order 
ones. In both cases the convergences of updating the mass function are rapid well for 
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Figure 9: The functional forms of the real and the imaginary parts of the 

mass function for \i along T = line. 



solving the Schwinger-Dyson equation. Unfortunately, near the phase transition line in the 
middle region, the convergences are too bad to obtain solutions with a suitable accuracy. 
However, an natural guess will be that the order of phase transitions continuously changes 
from first order to second order, through weak first order, in the middle region shown as 
in Fig. |10[ This type of diagram is also obtained in Refs. |Tl|, which is the same as 
that of two dimensional Gross-Neveu model in Refs. |21 




169 MeV 

Figure 10: The schematic view of the phase diagram from our result. 



4 Summary and Discussion 

In this paper we study the chiral symmetry restoration at finite temperature and chemical 
potential in QCD. We use the improved ladder approximation and the imaginary time 
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formalism. The improved ladder approximation does not violate the chiral symmetry, 
since the high energy behavior of the quark mass function is consistent with the result 
of the operator product expansion||19|| even at finite temperature and chemical potential. 
The phase transition point (or line) is determined by using the three order parameters; 
i.e., the VEV (if)ip) icev renormalized at 1 GeV, the quark mass gap S (0) and the pion 
decay constant f„. In the improved ladder approximation the infrared regularization pa- 
rameter tp must be introduced as in Eq. (|2.6| ) in order to regularize the running coupling. 
We, however, observe that the physical quantities do not depend on the parameter If- 
Then, our results are obtained without any degrees of freedom other than Aqcd which is 
determined by putting f n = 93 MeV at (T, fi) = (0, 0). 

In the case of the vanishing chemical potential /i = we have a second order phase 
transition at T c = 169 MeV. The critical exponents are extracted in Eq. ( |3.3|) . This 
shows that the QCD in the improved ladder approximation is different from mean field 
theories. In the case of the vanishing temperature T = we have a strong first order 
phase transition at \i = 598 MeV. 

In the (T, /i) = (0, 0) limit the functional forms A(p) = 1 and B(p) = as in Eq. ( |2.8|) 
give the solution of the Schwinger- Dyson equation in Landau gauge. In this paper we put 
these forms for any (T, /x) for simplicity. We should check the validity of this prescription. 
In the middle region (0 < T < T c and < fx < yU c ) near the phase transition line the 
calculation of the mass function is too hard for the error to vanish in the iteration method. 
It is necessary to obtain more efficient method for solving the Schwinger-Dyson equation 
in this region. After these problems being settled, the framework of the improved ladder 
approximation becomes a more convenient tool to figure out the nature of chiral symmetry, 
since it is easy to introduce fermions in the chiral limit. 
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